A RECURSION FOR DIVISOR FUNCTION OVER 
DIVISORS BELONGING TO A PRESCRIBED FINITE 
SEQUENCE OF POSITIVE INTEGERS AND A SOLUTION 
OF THE LAHIRI PROBLEM FOR DIVISOR FUNCTION 
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C") , Abstract. For a finite sequence of positive integers A = {a,}^ =1 , we 
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prove a recursion for divisor function ai (n) = J2d\n deA^ ■ As a 
corollary, we give an affirmative solution of the problem posed in 1969 
by D. B. Lahiri [3]: to find an identity for divisor function cr x (n) similar 
to the classic pentagonal identity in case of x = 1. 
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1. Introduction and main results 
We start with the two well known beautiful classical recursions. Let pin) 

m : 

be the number of all partitions of positive integer n and er(n) be the sum of 
its divisors. Then (sf [1],[5]) we have 

(1) p(n) = p(n — l)+p(n — 2)—p(n — 5)—p(n — 7)+p(n — 12)+p(n — 15) — ... 

(2) a(n) = a(n-l)+a(n-2)-a(n-5)-a(n-7)+a(n-12)+a(n-15)-... 



^ ■ where the numbers 1,2,5,7,12,15,... appearing in the successive terms in 

c3 ' are the positive pentagonal numbers {v m } given by 

(3) v m = m(3m =)= l)/2, to =1,2,... 

In identities (l)-(2) we accept that p(m) = 0, a{m) = when to < 0. The 
only formal difference is that (1) is true with the understanding that 

(4) p(0) = 1, 
while (2) is valid with the understanding that 

(5) o-(O) = n. 

Note that, formulas (l)-(2) are proved with help of the famous Euler pen- 
tagonal identity 

oo oo 

(6) JJ(l- g »)= {-l) m q m(3m - 1)/2 ■ 

71=1 771= — OO 
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In its turn, a combinatorial proof of (6) is based on the following statement 
(sf [1]). Let p e (n) (p (n)) denote the number of partitions of n into even 
(odd) number of distinct parts. Then 



Let <y x {n) denote the sum of the xth powers of the divisors of n. In 1969, 
Lahiri [3] noticed that every definition of Ofc(0) = f(n), k ^ 1 is irrelevant 
in order to keep the classical identity (2) and posed the following problem: 
"Whether analogous identities exist for divisor function Ofc(n) of higher 
degree?" Formally, for every not necessarily integer value of x, —oo<x< 
oo, for cr x (n) we could consider an identity of type (2) of the form 

a{n) = g x (n) + a{n — 1) + <r(n — 2) — a{n — 5) — a{n — 7) + 

where {g x {n)} is some "compensating sequence," and a solution of the Lahiri 
problem consists of a description of the compensating sequence for every n 
without a reference to its divisors. In particular, by the definition of cr x (n), 
and accepting as in (5) cr(0) = n, we find 



g x (l) = 0, g x (2) =r-2, g x (3) = 3* - 2* - 1, g x (4) = 4* - 3* - 1, 
g x (5) = 5 X -A X -3 X -2 X + A, g x (6) = 6 X - 5 X - A x , ... 



At first sight, this sequence is even more complicated than cr x (n), and it 
seems hardly probable to find a required description of it. Our paper, in 
particular, is devoted to this aim. For a simplification of our transforma- 
tions, below we accept the unique convention 



It is easy to see that in this case we have only a little change of the com- 
pensating sequence in the identity of the same form 

a x (n) = h x (n) + a x (n - 1) + a x (n - 2) - a x (n - 5) - a x (n - 7) + 

such that 



Note that this relation is so simple only due to Euler pentagonal identity 
(6); in more general case (see below Theorem 1) the corresponding relations 
could be very complicated and the convention (8) plays the unique role for 
the obtaining of general result. In particular, we write (l)-(2) in just a little 
another form. Namely, according to (8), instead of conventions (4)- (5), we 
accept the unique convention 





(8) 



a x {n) =0, if n < 




p(0) = 0, (7(0) = 0. 
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Then with help of (7) it is easy to see that, instead of (l)-(2), we have 

p(n) = h ip \n) + 

(9) p(n-l)+ p(n - 2) - p(n - 5) - p(n - 7) + p(n - 12) +p(n - 15) - 
where the compensating sequence has the form 

(10) = I'- 1 '™" 1 ' «/» = ™(*»=Fl)A . 

I 0, otherwise 

and, in view of the same structure of (1) and (2) and taking into account 
(4)-(5), we see that 

a(n) = h i - <T \n) + 

(11) a(n - 1) + a(n - 2) - <r(n - 5) - <r(n - 7) + a(n - 12) + a(n - 15) - 
where 

-l)" 1 - 1 ^, i/ ra = m(3m =F l)/2, 
0, otherwise 

Before formulating a generalization of (9) and (11), we study the divisor 
function over divisors belonging to a prescribed finite sequence A of positive 
integers. In the trivial case of a one-element sequence A = {a} we put 

(13) = ff' " >0, • * 6 (-«,,+«,). 

10, otherwise 

According to (13), we accept 

(14) <y [ x a] \n) =0, n < 0, 
such that 

(15) a«-> ) (n) = ^ a}) (n-a) + (f *{ " = °' . 

I (J, otherwise 

Consider now, for a fixed /c > 1, an arbitrary sequence 



(12) /i (<T) (n) 



(16) A = {a,-}} =1 

of positive integers. For a fixed x, let us consider an associated sequence 

(17) B(A;x) = {b l (x)}f =1 , 
where 

(18) bi(x) = a x h + a x n + a* + ... + a* , 
if the binary expansion of i — 1 is 

(19) i - 1 = 2 jW + 2 jW + ... + 2 jW , 1 < ji < j 2 < ... < j r , 1 < r < fc. 
In particular, since 2 fc — 1 = 2 1 " 1 + 2 2_1 + ... + 2 fe_1 , then 
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(20) b 2 k(x) = a x + a x + ... + a x k , 

while, since to i — 1 corresponds the empty set of terms in (19), then 

(21) b l (x) = 0. 
Furthermore, 

(22) b 2 (x) = a x , b 3 (x) = a 2 , h(x) = a\ + a 2 , etc. 
Moreover, denote 

(23) bi(l) = bi, l<i<2 k . 
For n > 1, consider divisor function over sequence A 

(24) *i A) (n)= d * 

d\n, deA 

in the understanding that every term d x repeats correspondingly to the 
multiplicity of d in sequence A. Besides, we accept the convention 

(25) <?i A \n) = 0, if n < 0. 

Denote by {£„} the Thue- Morse sequence [4], [2] which is defined as 

(26) t n = (-l)'W 

where s(n) denotes the number of ones in the binary expansion of n. 
Theorem 1. In convention a(n < 0) = 0, we have the following recursion 

(27) *i A \n) = hi A \n) + ^t 2l ^\n - h) 

where the compensating sequence h^\n) is defined as 

(28) hi A \n)= J2 tn-Mx)- 

i>2: bi=n 

Remark 1. Taking into account that 

1 + s (i - 1) = s(2(i - 1) + 1) = s(2i - 1), 
we prefer to write t 2 «-i instead of — 

Note that, as follows from (28), for n> b 2 k, hi A \n) = such that 

2 k 

(29) ai A \n) = Y^t 2l ^ A \n - b t ), n > b 2k . 
Consider now the divisor function 
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(30) a x (n) = J2d X - 

d\n 

Putting here 

(31) bi(x) = ft + f 2 + 3l + ... + f r (and b t = b t (l)) 

if the binary expansion of i — 1 is defined by (19), we obtain the following 
result. 

Theorem 2. We have 

<j x (n) = h x (n) + 

(32) a x (n-l)+a x (n-2)-a x (n-5)-a x (n-7)+a x (n-12)+a x (n-15)-..., 
where the compensating sequence {h x (n)} is defined as 

(33) h x (n) = ^2 t 2 i-ibi(x), n>\. 

i>2: bi=n 

Theorem 2 gives a solution of the Lahiri problem for divisor function 
a x (n). 

2. Proof of Theorem 1 

We use the induction over the number of elements of sequence A, the 
base of which is given by (15). Note that if, instead of A = {a 1; ...,ak}, to 
consider the sequence 

(34) A' = {ai, a k , a k +i}, 
then we have 

(35) ^ (n)=^W + ^ l}) (n). 

Furthermore, in the case of A', to every i, 1 < i < 2 k , with the binary 
expansion (19) of i — 1 corresponds bijectively the number 2 k + i from 
[2 k + 1, 2 k+1 ] with the expansion 

2 k + i - 1 = 2 jW + 2 jW + ... + 2 jr ~ 1 + 2 k 

such that the associated sequence has the form 

< + d% + a% + ... + a* , if 1 < i < 2 k , 



(36) h(x) 



a% + a% + al + ... + a* + a% +1 , if 2 k + 1 < i < 2 k+1 . 



This means that, for 1 < / < 2 k , we have 

(37) bi +2 k(x) = bi(x) + al +1 (in particular, b 1+2 k(x) = a£ +1 ). 
Notice also, that 
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(38) t 2 k+i +1 — 1; t2(z+2 fe )-i — ^2«+2 fc + 1 -i — —hi-i 
and 

(39) *2i-i = - E = 0. 

l<Z<2 fc l<K2 fc 

Suppose now that the theorem is true up to k. Then, using (37)- (38), we 
have 

2 fe+i 

i=2 

i=2 i=2 fc +l 

J^^i-i^'^n - bi) + 5^* 2 (i+2*)-i^- A ' ) (n - 6 i+2 fc) = 

i=2 i=l 
i=2 1=1 

Furthermore, by (40) and (35), we have 

2 fc+i 

(n-&0 = 

i=2 

2 fc 2 fc 

J>*-i^(™ - bi) + $>,-i<x« afc + l}) (™ - bi) 

i=2 i=2 
2 k 2 k 

(41) - J2t 2l ^ai A \(n - b t ) - a k+1 ) - 6,) - a fc+1 ). 

i=l i=l 

Note that, according to (15), 

2 fe 2 fc 

5>,-i<4 K+l}) (™ - 6,) - J]t2i-i<7« 0fc+l}) ((« " h) - a fc+1 ) 

i=2 i=l 

(42) =4 {ak+l}) (n) + at +1 £ i^. 

l<i<2 fc : n— bi=afe + i 

Therefore, from (41) we find 
2 k+i 

»=2 l<j<2 fe : ra-6i=a fc+ i 
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2 k 2 k 

(43) ^2t 2i -i<Ti A \n - h) -J2t 2l _ l( ri A \(n - a k+1 ) - b t ), 

i=2 i=l 

or, using the inductive hypothesis, we have 

2 fc+i 

i=2 l<i<2 fe : n-b 4 =a fe+ i 

(44) {a x A \{n - a k+1 ) - h {A \n - a k+1 )) + a x A \n) - h^ A \n). 
Furthermore, 

E t 2i -ibi(x) = E t 2i -ibi(x)+ 

2<i<2 k + 1 : bi=n 2<i<2 k : bi=n 

E hi-iHx) = h (A) - E *2i-i&j+2*(a0 = 

2 fc +l<i<2'=+ 1 : fe i= n l<Z<2 fc : fe, +2 fc =n 

(45) ^) (n) _ /l (A) (n _ afe+i) _ <+i £ ^ 

l<Z<2 fe : bi=n—ak+i 

Finally, summing the results of (44) and (45), we complete our proof: 

2 k+i 

^t 2l _ l( j {A '\n -h)+ E hi-Mx) = 

i=2 2<i<2 k + 1 : bi=n 

a x ^\n)+a x A \n) = a x A '\n)M 
3. Proof of Theorem 2 

If to consider as a finite sequence A the sequence A — A k — {1, 2, k}, 
then, for n < k, we have 

(46) 4 Afc) W = 

and, by Theorem 1, the (±)-structure of a x Ak \n) is the same as in the case 
of x — 1 (see (11)). Therefore, independently from the summands (either 
(Ji{n) or (7 x (n)) we have the same reductions, i.e. 

a x (n) = ai A *\n) = h^\n) + 

a x Ak \n - 1) + ^(n - 2) - <7< A *>(n - 5) - a x Ak \n - 7) + 

(47) a x Ak) (n - 12) + a x Ak \n- 15) - (n<k), 
with the compensating sequence 

(48) h x A *\n)= £ ^.^(x) 

j>2: fei=n 
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where bi(x) are defined by (31). If, instead of Ak, to consider N, then 
for every n we actually consider a finite part of (47) which corresponds to 
A n = {1,2, ...,n}. Thus (47) is true for A = N, and (32)-(33) follow.B 

Example 1. Consider the case of x = 1, i.e. the case of sum- of- divisors 
function. 

Then we have 

h[ N \n)=n J2 *«-i = -n £ (-l)*" 1 ' = n ( Po (n) - p e (n)) 

i>2: bi=n i>2: bi=n 

and, in view of (7), we obtain (11) as a special case of Theorem 2. 

4. Expression of compensating sequence {h x (n)} = {h x (nY N ^} via 

KNOWN SEQUENCES 

Note that from the definition of sequence b n (x) (see (31) and (19)) it 
follows that if 

(49) n-l = ^(i)*" 1 

i>l 

is the binary expansion of n — 1, then 

(50) bn(x) = Y,j3(i)i x , 

i>l 

such that 

(51) 6 n = J]/3(i)i. 

i>l 

Notice that, (51) is Sequence A029931(n-1) in [6]). Denoting 

(52) A029931(n) = r/(n), 
according to (33), we have 

(53) h x (n)= (-l)^-\ +1 (x). 

Example 2. Consider the case of x = 0, i.e. £/ie case of the number of 
divisors of n. 

Then, by (50) and (53), the compensating sequence has the form 

(54) ho(n)= J2 (-lY^-'sU), 

J>1: V(j)= n 

where s(n), as in the above, is the number of ones in the binary expansion 
of n. The first terms of compensating sequence {h^\n)} n >i are: 
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(55) 1, 1, -1, -1, -3, 0, -2, 1, 2, 1, 2, 4, 1, -1, ... 

Note that, in view of (49), the expression (51) gives the number of all 
partitions with distinct parts of a fixed values of b n . This means that if to 
denote by e n the set of the terms of A029931 for which rj(j) = n : 

ei = {1}, e 2 = {2}, e 3 = {3, 4}, e 4 = {5, 8}, e 5 = {6, 9, 16}, 
e 6 = {7, 10, 17, 32}, e 7 = {11, 12, 18.33.64} 
then the concatenation of this sets leads to the ordering of all partitions of 
n with distinct parts : 

n = ^2(3(i)i 

i>l 

respectively to the values of ^2i i>1 f3{i)2 l ~ 1 . Thus this way leads us to the 
Adams- Watters sequence " Decimal equivalent of binary encoding of parti- 
tions into distinct parts" (see A118462 in [6]): 

(56) 1, 2, 3, 4, 5, 8, 6, 9, 16, 7, 10, 17, 32, 11, 12, 18, 33, 64, ... 

Denote Sequence (56) via W(n) and put |e n | = R(n), where, for n > 
1, {R(n)} is Sequence A000009[6] that is the number of partitions of n 
into distinct parts. Finally, denote 

n 

(57) T(n) = ^2R(k) = A036469(n) - 1, n > 1. 

k=i 

Then from (53) we find 

R(n)-1 

(58) h x (n) = J2 (-l) s(H/(T(ri) " m)Kl ^(T(n)- m)+ i(^)- 

m=0 

Example 3. Let us calculate the seventh term h (7) of sequence (55). 
By (58), here we have 

R(n)-1 

(59) ho(n)= ^ {-l) s{w{T{n) - m)) - l s{W{T{n)-m)). 

m=0 

If n = 7, then we have from the corresponding tables of [6]: 

R(n) = 5, T{n) = 18, W(18) = 64, W(17) = 33, 
W(1Q) = 18, W(15) = 12, W(U) = 11. 
Thus, according to (59), we find 

h (7) = 1- 2-2-2 + 3 = -2. 
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5. Some another identities 
In case of the finite set 

A k = {2^} k j=1 , 

according to (18), we have 

(60) bi(x) = 2 {jl ~ l)x + 2 {j2 ~ 1)x + ... + 2^" 1)x , 
if 

(61) i-l = 2 jW + 2- ?w + ... + 2 ir ~\ l<ji< j 2 < ■■■ <j r , 1 < r < k. 
Thus bi has an especially simple form: 

(62) bi = i-l. 
Let 

n = 2 ai ~ l + ... + 2 am ~ 1 . 
Then, according to Theorem 1 and (62), we find 

h { x Ak \n)= t2i-Mx) = 

i>2: bi=n 

(63) *2n+A+l(aO = (_l)-(2n+l)( 2 K-D- + ... + 2 ^ l >) 

and 

a i Ak \ n ) — (-l) s(2ri+1) (2 (ai ~ 1):E + ... + 2 ( - am ~ 1 ^ x ) + 

(64) £ ( _i)<<*-i)^) (n _(<_!)). 

i>2 

Considering now the infinite sequence of powers of 2: 

A = {2 J - 1 L> 1 , 

we conclude that 

a { x A \n) = (_i)«(2n+i)( 2 ( Q? i- 1 ^ + ... + 2 {am ^ x ) + 

n 



i>2 

or 



al A \n) = (— l) s ( 2n+1 ) {2^ ai ^ x + ... + 2 ( -° lm ^ x ) — 

n-l 

(65) E(- 1 ) S(n "' ) ^ ) (i)), 

j'=i 

where 

n = 2 ai ~ l + ... + 2 am ~ 1 . 
In particular, in the case of x — 0, we obtain the identity 

n 

(66) E(- 1 ) S( "" j) ^ A) 0')) = (-l)^" 1 *^). 
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The sequence {<Jq (n)) — l} n >i is well-known so-called "the binary carry 
sequence" (A007814 in [6]). In the case of x — 1, we obtain the identity 



The sequence {<r\ (n))} n >i is also well-known (see A038712 in [6]). 
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